Geometric Phase in PT-Symmetric Quantum Mechanics 
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Unitary evolution in 'PT-syninietric quantum mechanics with a time-dependent metric is found 
to yield a new class of adiabatic processes. As an explicit example, a Berry-like phase associated 
with a "PT-symmetric two-level system is derived and interpreted as the flux of a fictitious monopole 
with a tunable charge plus a singular string component with non-trivial phase contributions. The 
Hermitian analog of our results is also discussed. 
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Parity-time-reversal (VT) symmetric quantum me- 
chanics (QM) has emerged as a complex extension of 
conventional QM, via non- Hermitian Hamiltonians that 
possess unbroken PT-symmetry and hence real energy 
spectrum On a deep level T'T-symmctric QM is 

not expected to challenge the current quantum theory. 
Nonetheless, some fascinating features of PT-symmetric 
QM have been revealed , with their implications being 
pursued both experimentally Q and theoretically It 
is now certain that PT-symmetric QM will bring more 
insights and results that are not obvious in conventional 
QM. 

In PT-symmetric QM, a given non-Hermitian Hamil- 
tonian iJ(X), here parameterized by a multi-dimensional 
vector X, is sclf-adjoint with respect to a generalized in- 
ner product. In particular, a Hermitian, positive-definite 
metric operator WCK) constructed from the VT sym- 
metry of a system defines a new inner product, denoted 
(^', for two arbitrary states ^' and <f>. In Dirac's 
bra and ket notation, (*, = (*|M^(X)|$). The sclf- 
adjoint condition of H(X) then becomes 



W{X)H{X) = H^{X)WiX 



(1) 



which yields real spectrum and orthonormal eigenstates. 

Because the metric operator W{X) must be always at- 
tached to H{X), a time-dependent H{X) via X = X{t) 
calls for a time-dependent metric, an interesting situa- 
tion previously overlooked. This motivates us to exploit 
time-dependent problems in PT-symmetric QM to ex- 
plore fundamental issues, among which adiabatic manip- 
ulation and Berry-phase related problems are of broad 
interest [1, Indeed, to preserve the W(X)-based inner 
product and hence unitarity, the time-evolution equation 
in PT-symmetric cases is shown to have special features, 
thereby defining new dynamical problems from the very 
beginning. 

Specifically, we present below (i) general considerations 
of the adiabatic time evolution associated with a time- 
dependent H[X{t)] in PT-symmetric QM, and (ii) in- 
triguing geometric phase results based on a six-parameter 



complete characterization of 2 x 2 PT-symmetric Hamil- 
tonian matrices ■ The found Berry-like phase in a rep- 
resentative case is interpreted as the flux of a fictitious 
monopole with a tunable charge, plus a singular com- 
ponent similar to a Dirac string but with a non-trivial 
effect. To gain more insights, we also analyze the Hermi- 
tian analog of our non-Hermitian problem, finding that 
the context of time-dependent PT-symmctric QM also 
opens up an interesting class of geometric-phase prob- 
lems in conventional QM. Furthermore, we discuss com- 
mon features shared by this work and a recent study 
regarding geometric phases in nonlinear systems. 

Consider then the following linear evolution equation 
for PT-symmetric QM (setting h ~ 1 throughout). 



(2) 



The general form of A(t) can be determined by the unitar- 
ity condition. That is, for two arbitrary solutions |$(i)) 
and !«'(*)), 



dt 



{m\w[xitm{t)) = o. 



(3) 



Equations ^ and (jS]) directly lead to 
.dW[X{t)] 



dt 



= AHt)W[Xit)]-W[Xit)]A{t). (4) 



Letting 



A{t) = H[X{t)] - tK[X{t)], 



(5) 



where H[X{t)] is the time-dependent PT-symmetric 
Hamiltonian satisfying Eq. ([T]), one arrives at the fol- 
lowing condition for K[X{t)]: 

VW[Xit)] ■ ^ = K^[Xit)]W[Xit)] + W[Xit)]K[Xit)], 

(6) 

where V is the gradient in the X-space. Several re- 
marks on Eq. ^ are in order. First, if K[X{t)] satisfies 
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Eq. del, so does K[X{t)] + iH{t), where H{t) is an arbi- 
trary T^T-symmetric Hamiltonian with W['X.{t)]H{t) ~ 
H''' {t)W[K{t)]. This arbitrariness hence just reflects a 
different partition of A{t). Second, in general K[X.{t)] is 
a function of cDL/dt, the velocity in the X-space. Third, 
to find an interesting sample solution to Eq. ([5]), we 
put forward below an additional condition K[lL(t)] = 
i4r^[X(<)]. With these clarified, one can now adopt a so- 
lution K[X{t)] = M[X{t)] ■ (dX/dt), with 

VW{X) = M(X)M^(X) + T^(X)M(X) (7) 

and M(X) = Mt(X). Clearly, if dX/dt 0, then the 
constructed K[X{t)] — > 0. 

Next we examine the time dependence of [see 
Eq. ([2])] for a slowly-varying X. Because in the adia- 
batic limit, A{t) = H[X{t)], the natural adiabatic rep- 
resentation should comprise instantaneous eigenstates of 
H[X{t)], denoted |?/'„[X(<)]). Specifically, if 



H[X{t)] |Vn[X(<)]) = E,,[X{t)] |V„[X(t)]) 



(8) 



with the orthonormal condition 

(V'„[X(0]|W^[X(i)]|V'™[X(t)]) = Srnn, then we can 
expand |\E'(t)) as 



(9) 



where c„ is the expansion coefficient, and a„ is a dy- 
namical phase defined as the time-integral of the instan- 
taneous real energy eigenvalue En[X{t)], i.e., an{t) = 
- J* En[X{t')]dt'. Substituting Eq. (O into Eq. ([2]), one 
obtains 

where 



Trnn (^) 



(V™(X)|W^(X)|VV™(X)) 



(V'™(X)|W^(X)M(X)|^™(X)) } ■ — , (11) 



dX 



(V™(X)|W^(X) 



VH{X) 



EniX)-E,niX) 



l^n(X)) 



dX 



+ (V'™(X)|H^(X)M(X)]|V4X))| • ^ (12) 

with X = X{t). Clearly, due to the oscil- 
lating phase factor e*["" in Eq. (UHl), the 
transition probabilities between non-degenerate eigen- 
states |-0„(X)) and |'(/'„i(X)) are negligible if we have 
VH(X) ■ {dX/dt)/[En{X) - E^{X)Y < 1 and M(X) • 
{dX/dt)/[En{X) — ETn{X)\ ^ 1 (a precise adiabatic con- 
dition is out of the scope of this work and the analo- 
gous topic in conventional QM is still under study). As 



such, in the limit of dX/dt — )• 0, adiabatic evolution in 
T'T-symmetric QM does exist even for a time-dependent 
metric, insofar as the time evolving state can remain to 
be one instantaneous eigcnstate of a time-evolving Hamil- 
tonian. 

For adiabatic evolution, the Qm term in Eq. pUj) gives 
the solution c,„(t) = Cm(0)c'''''*-', with 

lit) - * r'\{4>ux)\w{x)\v^^{x)) 

Jx(o) 

+ (V™(X)|M/(X)M(X)|Va„(X))] • dX. (13) 

Because ^{t) is solely determined by the geometry of a 
navigation path in the X-spacc and does not depend on 
the duration of the adiabatic process, it is of a geometric 
origin. In particular, for a cyclic path, this geometric 
phase becomes a Berry-like phase 7^, i-e.. 



IB ^ if [(7/'™(X)|M^(X)|VV;™(X)) 

(^„(X)|W^(X)M(X)|V'„(X))].dX. (14) 



In Eq. p4|) , the first term is parallel to a familiar Berry- 
connection integral 0, but the second term represents 
effects due to a time-dependent metric. It is the sum 
of these two terms that can ensure a real 75. Note 
also that 75 is invariant upon the gauge transformation 
\ipm{X)) e*^(^)|V'm(X)), where /(X) is an arbitrary 
continuous phase function. 

Let us now apply these general considerations to a two- 
dimensional Hilbert space. Here, the T operator is taken 
as the Dirac conjugate, and a nontrivial 2x27^ operator 
is chosen as 



cos f 
sin6' c" 



sni( 



cos ( 



(15) 



where a = (cr^;, ay, CTz) represents the standard Pauli 
matrices, and eV = (sin cos sin0 sin <p, cos 0). To 
avoid confusion with the vectors in the X-space of ar- 
bitrary dimension, below we always use A to represent 
a three-dimensional vector A. The general form of VT- 
symmetric 2x2 Hamiltonian matrices can then be spec- 
ified by six real parameters [X = (e, a, 6, 9, 0, 6)] 0: 



H 



2x2(X) = e + [aer + ibcosSee + ibsinSe^p) ■ a, (16) 



where eg = (cos0 cos (p, cos sin (/s, — sin0) and = 
(— sint/?, cos 0). Here we apply the constraint > 6^ 
to ensure a non-degenerate real spectrum (see explicit 
eigenvalue expressions below). The special case of 6 = 
recovers Hermitian matrices with four real parameters. 
Attached to such a complete characterization of i/2x2(X) 
is the following 2x2 metric operator, 



T^2x2(X) 



V^2— 52 



l + /3-a 



(17) 
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where 



P = — (sin 5 €0 — cos 5e^p). 



(18) 



The real eigenvalues of H2X2OQ arc E± = e± \^a'^ — b'^, 
with eigenvectors 



1 

A4 



e "f" (a sin 6' + z6 cos (5 cos 61 + 6 sin (S) 
—a cos 9 + ib cos i5 sin 9 ± \/ a? — b"^ 



(19) 



Here the normalization factor A/± is determined by M\ 
2\a\\J d?- — h"^ [1 + (6/a) sinJsinfl =F W cos0], where 



- 52 



(20) 



The fact that i72x2(X) has two more free parameters 
than a 2 X 2 Hermitian matrix is stimulating. In the 
following we focus on those navigation paths that only 
involve changes in 9 and (j). Extension to other cases (e.g., 
change in 5) is certainly worthwhile but straightforward. 
In the (0, 0) manifold, a, 6, and \ j3\ arc constants and we 
obtain a simple Hermitian solution to Eq. ([7]): 



M2x2(X) = iV^(X).a. 



(21) 



Substituting Eqs. (HZ]), ([19]), and ([ST]) into Eq. dH]) leads 
us to a lengthy calculation, but we finally arrive at a 
rather compact expression for the geometric phase. For 
example, associated with the Bcrry-likc phase de- 

noted 7^ is found to be 



-f+ = [F^{X)dcl, + Fe{X)de] , 



where 



{b/a) cos 5 



2[1 + {b/a) sin SsmO-U cos 9] ' 



(22) 

(23) 
(24) 



To better understand 7^ we now treat 9 and as two 
spherical angular coordinates on a sphere defined by r = 
rcr = (a;, y, z). If we define 



A 



rsm9 



(25) 



then because F^{X)d(l) + Fe{X)d9 
becomes 



A ■ dr. Eq. 



it 



A ■ dr. 



(26) 



This motivates us to regard A defined above as a three- 
dimensional vector potential in the r-spacc. If the sur- 
face enclosed by a closed path on the sphere does not 
include = 0,7r, i.e., the north or south pole, then A is 



well behaved everywhere on this surface, and as a result 
Stokes' theorem converts Eq. ((26)) into a surface integral: 
1b ^ J ^ A - dS. Further using dFg{X)/d(j) = 
and the standard curl operator in spherical coordinates, 
one finds 7^ = -{U/2) ] J sm9d9d(l) = -{U/2)n, where 
il is the solid angle spanned by the closed path. On the 
other hand, if a closed path does enclose the north pole of 
the sphere (for convenience we always assume a counter- 
clockwise path) , then the circulation integral in Eq. ([25)1 
can be converted into two integrals: one is on a deformed 
closed path circumventing the north pole and the other is 
on an infinitely small circle around the north pole. In this 
second situation, we obtain 7^ = —{U/2)n + (1 +U)tt. 
The sole contribution by the north (or south) pole can 
then be found by setting 51 = (or ft ~ An) . This is con- 
sistent with a direct calculation from Eq. (j22[) by fixing 
6* at (or tt). 

With a fictitious magnetic field denoted Bcs, for 
an arbitrary closed path on the sphere can now be ex- 
pressed as a magnetic flux through a surface enclosed by 
the closed path, i.e., 7^ = / / Bcff • dS, where 



4ff = {l+U)Tr6{x)d{y)e, 



U r 
2^ 



(27) 



with 6*2 = er{9 ~ Q) . Equation (|27p is our main result for 
a Berry-like phase of a T^T-symmetric two-level system. 
-Beff represents a fictitious singular field pointing at the 
north pole, plus a virtual magnetic monopole. The mag- 
nitude of the charge of such a virtual monopole is con- 
tinuously tunable: it is given by a unit charge go = ^\ 
multiplied by the factor U £ (0, 1] for a > or e [-1,0) 
for a < 0. The singular component of Eq. (P7)) is anal- 
ogous to the so-called Dirac string in the Dirac quanti- 
zation rule for the product of a magnetic charge and an 
electric charge. Remarkably, the contribution by the sin- 
gular string component here is seen to be physically rel- 
evant: it generates a flux in the range of [0, tt) or (tt, 2tt]. 
Note that, in the limit of 6 = and hence U = ±1, the 
monopole fleld in Eq. (P?]) carries a charge igo- Under 
the same limit the flrst component of Eq. (|27p reduces 
to the trivial and unobservable Dirac string because it 
can only contribute a phase of or 27r. This makes it 
clear that the standard results for the Berry phase of a 
Hermitian two- level system can be all recovered ioi b ~ 0. 

The issue of geometric phase in non-Hermitian sys- 
tems was first touched upon in a very interesting study 
in Ref. However, the spectrum studied therein is not 
real. In addition, the construction of a geometric phase 
in Ref. Q was based upon a time-dependent Schrodinger 
equation and its non-equivalent adjoint. By contrast, 
here the spectrum is always real, the adiabatic condition 
is similar to that in the Hermitian case, and the evolu- 
tion of any state in the Hilbcrt space satisfies a common 
dynamical equation. 

To shed more light on our results we next examine a 
Hermitian analog of our geometric-phase problem. In- 
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deed, the real spectrum of a PT-symmetric Hamiltonian 
^^2x2(X) implies the existence of a Hermitian Hamilto- 
nian /i(X) with the same spectrum [§]. In particular, 
/i(X) can be constructed from by the following 

similarity transformation. 



MX) =77(X)F2X2(X)77-1(X), 



(28) 



where the transformation matrix 77(X) satisfies 7?^(X 
W^2x2(X). One explicit form of 77 (X) is found to be 



2{\a\ + V^^^) 



[VF2x2(X) + l]. (29) 



Equations and lead to h{X) = e+U\a\er-a. It is 
noteworthy that /i(X) is seen to be a standard Hermitian 
Hamiltonian describing a spin-i particle in a magnetic 
field. 

Returning to Eq. ^ using the representation of 
\ipr,{t)) = rj(X.)\^>{t)), one finds 

= hm)]\A,it))+v[X{t)]\4'4t)), (30) 

where the Hermitian potential function w(X) is given by 
«(X) = irjiX) [rj-\X)VviX) 



-M2x2(X)] rri(X) 



dX 



(31) 



Note that if M2x2(X) = 7j-\X)Vt]{X), then t;(X) = 0. 
In this special case the dynamical equation for \ipr,(t)) 
is solely determined by the familiar Hamiltonian ft,(X). 
This fact further confirms that the dynamics associated 
with a time-varying metric in PT-symmetric QM is a 
legitimate extension of the stationary-metric cases con- 
sidered in Refs. [13, El- 

Of more significance, because M2x2(X) can be cho- 
sen to be Hermitian [see Eq. (|2T|) ] and r/^^(X) Vr/(X) 



is not Hermitian in general, we have 77^^(X)Vr/(X) 
M2x2(X) 7^ 0. As such, for each choice of M2x2(X) sat- 
isfying Eq. ([7]), we have a different realization of w(X) 
from Eq. (|3ip . An interesting class of time-dependent 
problems in conventional QM can thus be defined from 
a PT-symmetric framework, always producing a Hermi- 
tian perturbative term dependent upon dX/dt, the veloc- 
ity in the parameter space. Within the conventional QM, 
such type of perturbation may also arise from a designed 
feedback mechanism. Alternatively, so long as a dynam- 
ical model (could be phenomenological, e.g., a system 
under a mixed quantum-classical description) requires a 
dX/dt-dependent effective Hamiltonian, then a term sim- 
ilar to v{X) emerges and all other terms in higher orders 
of dX/dt can be safely neglected in adiabatic processes. 

It is also evident why the nonzero v(X) term can be 
important for adiabatic processes. In the adiabatic limit, 
though the perturbation v{X) is vanishingly small, its ef- 
fect will be accumulated over a diverging time scale. In 



brief, the expectation value of v(X) on the instantaneous 
eigenstate, which is of the order of dX/dt, constantly 
corrects the quantum phase. The accumulation of this 
phase correction gives a nonzero and finite extra phase 
of a geometric nature. This further justifies the emer- 
gence of a non-conventional Berry-like phase as well as a 
non-trivial 7^ even at 9 = 0, tt. This "Hermitian" per- 
spective is also verified by our calculations based directly 
on Eq. ([30|). Indeed, due to the mapping from Eq. ([2]) 
to Eq. (|30l). an adiabatic process for |^'(i)) is mapped to 
that for \ijjr]{t)), and the Berry- like phase for |^(i)) is the 
same as that for \ipri{t))- 

In studies of the adiabatic evolution associated with a 



Gross-Pitaeviski equation, recently Fu and Liu [12| also 
found a remarkable Berry-like phase. In their model 
the instantaneous eigenstate changes adiabatically and 
at the same time induces a perturbation through a state- 
dependent mean-field Hamiltonian, with the perturba- 
tion strength considered to the first order of the velocity 
in the parameter space. This novel feedback mechanism 
clearly belongs to, and is also one explicit example of, 
the new class of adiabatic evolution found here. 

To conclude, adiabatic unitary evolution with a time- 
dependent metric is shown to exist in P7~-symmetric 
QM. For a T'T-symmetric two-level system, a Berry-like 
phase is found to display unusual features when com- 
pared with current notion of the magnetic monopolc and 
the Dirac string. Our results establish an entire class of 
new geometric phase problems for both T^T-symmctric 
QM and conventional QM, where adiabatic manipulation 
yields additional phase corrections to the conventional 
Berry-connection integral. This work therefore repre- 
sents another interesting application of 7^7~-symmetric 
QM with unexpected findings. J.G. is funded by NUS 
"YIA" (WBS grant No. R-144-000-195-101). Helpful 
discussions with Jie Liu, Li-Bin Fu, and Qi Zhang are 
gratefully acknowledged. 
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